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SUMMARY 


The inf luence of solid-body rotation on a screen-produced turbulence 
in a flow between concentric, rotating cylinders has been investigated. 
Radial distributions of the three components of turbulence intensity and 
the three turbulent shear stresses were measured at a fixed distance down- 
stream of the screen both with and without rotation. An energy balance 
is made for the three intensities from measured decay and dissipation 
terms. 

In addition, the change of the axial and peripheral mean velocity 
profiles with distance along the axis of rotation is presented. 


INTRODUCTION 


Although the bulk of theoretical work on turbulence has been 
restricted to homogeneous isotropic turbulence, there have appeared in 
recent years a number of investigations on turbulent shear flows char- 
acterized by the presence of mean velocity gradients. The ordinary types 
of shear flows that have been investigated, for example, channel (ref. l) , 
pipe (ref. 2), wake (ref. 3) , boundary layer (refs. 4 to 6), and jet 
(ref. 7), have one thing in common: The principal mean rate of strain is 

either exactly or approximately equal to the rotation, that is, one-half 
the vorticity, of the mean flow. A somewhat unusual experiment on a tur- 
bulent flow with uniform mean rate of strain but with no mean vorticity 
was recently reported by Townsend (ref. 8). The present investigation 
deals with the other extreme case, namely a turbulent flow with uniform 
mean vorticity but with no mean rate of strain. Such a flow was obtained 
by a solid-body rotation of the mean flow between two concentric circular 
walls, the walls also rotating like a rigid body about their common axis. 
An approximately uniform mean axial velocity served as a carrier. The 
turbulence was created by letting the flow pass through a series of 
screens, a customary method for generating turbulence except that in 
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this investigation the screens were attached rigidly to the tunnel walls 
and therefore also rotate with the walls. 

Even without the influence of viscosity and turbulence , the solid- 
body rotation of an ideal fluid has received considerable attention. 

(See a recent survey in ref. 9.) Taylor (ref. 10) showed that all small 
steady motions relative to the rotating fluid must be two-dimensional in 
the sense that any two fluid particles originally in a line parallel to 
the axis of rotation must remain so. This work was extended recently by 
Long (ref. 11). Morgan (ref. 12) considered the same problem without 
the restriction that the motion be steady. For the case of small forced 
oscillations that vary sinusoidally with time, he found that the differ- 
ential equations describing these oscillations are of elliptic type for 
cr > 2cd and of hyperbolic type for a < 2 a>, where a is the frequency 
of oscillation and od is the angular speed of rotation. He speculated 
therefore that, for a = 2co, the oscillation may grow to very large ampli- 
tudes. It was also found in the same reference that the ratio of ampli- 
tudes of the radial velocity to the tangential velocity is or/2o>; thus, 
for large values of oo compared with values of a, the radial component 
of the oscillations is suppressed, and the motion again appears to be 
two-dimensional. Stewartson (ref. 13) treated a weak source in a rotating 
fluid; the very peculiar result was obtained that the output of the source 
spreads axially but not radially, again a two-dimensional effect. 

An important feature of many flows encountered in practice (e.g., 
the flow in a turbomachine) is the fact that the streamlines may be 
curved, often quite strongly, introducing a pressure gradient in the 
direction perpendicular to the main flow direction. It is of great 
interest to know something of the effect that this curvature may have 
on the turbulence characteristics. One such investigation is reported 
in reference 14, where the effect of curvature on a fully developed 
channel flow is studied. It is clear that the curvature of the stream- 
lines alone is not sufficient to determine uniquely the characteristics 
of turbulence; the velocity distribution also must be considered. In 
the present investigation a solid-body rotation was therefore chosen for 
the simplicity of constant mean vorticity and zero laminar shearing 
stress . 

The present investigation was conducted at the Johns Hopkins University 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. The author gratefully acknowledges the 
important contributions and criticism of Dr. Hsuan Yeh, senior Faculty 
Advisor. He also wishes to thank Dr. Stanley Corrsin for his encourage- 
ment and suggestions. Dr. G. F. Wislicenus for his role in initiating 
the work, and Dr. Salamon Eskinazi, Dr. J. R. Weske, and Mr. W. G. Rose 
for numerous discussions throughout the investigation. 


NACA TN 4135 


3 



SYMBOIS 

A 

A A 

coefficient in relation — = 

■ § 

a 

radius of inner cylinder, 2.766 in. 

b 

radius of outer cylinder, 5-5^ in. 

c 

constant of integration 

d 

channel width 

E 

mean voltage across hot wire 

e 

fluctuating voltage across hot wire, dE 

F ^( K i) 

fraction of turbulent energy v^ associated with 
Ki, cm 3 /sec 2 

g. 

strength of disturbance defined in appendix B, v 0 /ax* 0 

i 

mean hot-wire heating current 

Io 

hot-wire heating current extrapolated to zero velocity 

K 1 

one -dimensional wave number in direction of mean flow, cnr-L 

L s 

integral scale of turbulence in direction of flow 

M 

mesh length, that is, distance between mesh wires (reciprocal 
of mesh number), in. 

N 

signal from hot wire defined in appendix A 

n 

coordinate measured perpendicular to streamwise direction on 
stream surface 

P 

mean static pressure 

P 

instantaneous pressure fluctuation 
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q. = 

R 

R o 

Re 

R xy 

R X 

r 

*i 

*0 


^v r 2 + v 0 2 + v z 2 

hot-wire resistance with heating 
hot-wire resistance with no heating 
Reynolds number, Vo(b - a)/v 
correlation coefficient, v x Vy/v x Vy 

turbulence Reynolds number, vA/v 

radial cylindrical coordinate 

radius of inner wall of plane curved channel 

radius of outer wall of plane curved channel 


S slope of hot-wire calibration curve 

s coordinate measured along streamwise direction 

U component of mean velocity parallel to surface 

V mean velocity component, direction indicated by subscripts r, 

0, z, n, and s 

V 0 mean throughflow velocity (V z averaged over flow area) 

Yp mean velocity component perpendicular to hot-wire axis 

v instantaneous velocity fluctuation, direction indicated by 

subscripts r, 9, z, s, n, x, and y 

v root-mean-square value of velocity fluctuation v 

X coefficient in the hot-wire response equation, 

sin p cos p sin a 

1 - sin 2 a sin^p 

X , y , z Cartesian coordinates 
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Y 


z 

a 

P 

6 

i 

n 

0 

A 

p 

v 

P 

a 

T l 


coefficient in the hot-wire response equation, 
sin a cos a sin 2 p 

1 - sin 2 a sin 2 p 

coordinate along axis of rotation; also axial distance from 
farthest downstream 20-mesh screen 

angle of turning of hot-wire probe, measured in plane per- 
pendicular to probe axis 

angle between hot-wire axis and radial direction 
boundary-layer thickness 

axial distance from farthest downstream 4-mesh screen 
exponent in relation I 2 = S^V s (cos a) 1 ! + I Q 2 

peripheral cylindrical coordinate 

2 v z 2 ? 2ve2 

microscale of turbulence, Av z , z = > A Vn . z = > 

and so forth 7~ To TZ To 

© © 

viscosity 

kinematic viscosity 
density 

frequency of forced oscillations 
laminar shearing stress 


T-j- turbulent shearing stress 

(f) flow angle , measured from the axial direction 

a) angular velocity of rotation 
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WIND TUNNEL 


The rotating wind tunnel was designed to produce a radially uniform 
axial velocity and solid-body peripheral velocity distribution, axially 
symmetric to prevent further complications. These combined requirements 
suggested the idea of using rotating screens. 

It is known that a flow, approaching with some angle to the normal 
to a screen, will be turned by the screen toward that normal so that 
downstream of the screen the velocity component parallel to the plane of 
the screen will be smaller than the corresponding component upstream of 
the screen. For a large number of screens, the parallel component of the 
approaching flow will be substantially reduced in the flow behind, this 
component vanishing in the limit of infinitely many screens . If the 
screens are rotated about a normal axis the flow behind them will also 
rotate like a solid body. Further, if the pressure-drop coefficient of 
the screens depends on the normal component of velocity only, as it does 
for approach angles of up to 45° (ref. 15) , nonuniformities in axial 
velocity will be diminished by the screens as well as peripheral 
nonuniformities . 


To avoid the excessively large number of screens that would be nec- 
essary to induce rotation by this method alone, an impeller with approxi- 
mately axial vanes at discharge was employed here to create a prerotation. 
This impeller was followed by a system of screens rigidly attached to it 
and thus rotating with it. These acted to reduce both nonuniformities 
in axial velocity and departures from solid-body rotation. They also 
acted to suppress peripheral variations caused by the wakes of the vanes. 


The impeller was a modified version of a previously existing mixed- 
flow runner (fig. l) . By machining, a vane system resulted whose outer 
surface was cylindrical and 10.688 inches in diameter, but whose hub 


varied in diameter from 



inches at the inlet to 5-532 inches at the 


outlet. Five screens, spaced l/4 inch apart, were attached to the hub 
in a single unit just downstream of the trailing edge of the vanes. 
They were 20-mesh screens, of 0.01- inch-diameter wire, resulting in a 
solidity of O.360 and a pressure drop coefficient of approximately 1.0 
for each screen for the throughflow velocity of 42.3 fps that has been 
used in these experiments. The screen Reynolds number based on rod 
diameter and throughflow velocity is 207 • As will be described, two 
additional turbulence -producing 4-mesh screens were later installed 


8 inches and 13^ inches, respectively, downstream from the last 20-mesh 

4 

The Reynolds number of these screens is 975* 


screen. 


NACA TN 4135 


7 


Since the impeller alone will not produce a sufficient throughflow, 
a 3. 75-horsepower axial blower is used to supply air to it at the rate 
of 4,000 cu ft/min through a settling chamber. A 20-mesh screen is 
located just before the impeller entrance. An adjustable bleed is pro- 
vided between blower and impeller to make up for blower characteristics 
unfavorable for this particular application. With the speeds of the 
blower and impeller separately controlled, variation in the ratio of 
peripheral to axial velocity is possible. 

Downstream of the impeller and screens, two coaxial cylinders form 
the inner and outer walls of the tunnel; the inner wall is rigidly 
attached to the hub of the impeller, the outer wall, rigidly attached 
to its shroud. 

The inner wall of the test section is formed by the outside surface 
of a piece of aluminum pipe 6 feet in length, with an outside diameter 
of 5*532 inches, supported by two bearings on a shaft extending through 
the inside of the pipe on which is also supported the impeller on its 
own two bearings. Both pipe and impeller rotate as a unit when the pipe 
is driven by a belt from the rear with a variable-speed direct -current 
motor . 

The outer wall of the test section consists of sections of various 
lengths of a cylindrical plastic tube having an inside diameter of 
10.688 inches and a wall thickness of l /2 inch, supported from the out- 
side on separate bearings. The ratio of inner to outer diameter of the 
test section is thus 5*532 inches to 10.688 inches or O. 518 . The distance 
between inner and outer walls is 2.578 inches. The length of this outer 
wall is adjustable by fitting together a different number of sections. 

As mentioned, the outer wall is rigidly attached to the impeller and 
thus also rotates with it. 

With impeller, screens, inner wall, and outer wall all rotating 
together as a unit, a special arrangement is necessary to introduce 
measuring probes into the flow. The scheme used is as follows: A 
stationary section of the outer wall is fitted to the downstream end of 
the tunnel. The measuring probe enters radially through a hole in this 
piece as close to the sealed joint between the rotating and stationary 
sections as possible. Thus the outer wall at the traversing station is 
not in rotation, but it is felt that the resulting sudden change in 
peripheral velocity is limited to a region very close to the wall. 

An assembly drawing of the wind tunnel is shown in figure 2. Fig- 
ure 3 shows in more detail the traversing section. 

The flow leaving the screens will undergo considerable change as it 
proceeds downstream. Thus, in contrast to investigations on fully devel- 
oped flow, several traverses are necessary at various distances from the 
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screens. This is accomplished by varying the length of the outer wall 
by fitting together the proper number of rotating sections and then 
adding the stationary traversing section at the end. 

Two general views of the tunnel are shown in figures 4 and 5> looking 
downstream and upstream, respectively. 


INSTRUMENTATION 


The instruments used for all mean flow measurements were a shielded 
total head probe (Kiel probe), a Pitot static tube, and a claw direction- 
finding probe. These instruments were mounted on a micrometer traversing 
apparatus to allow for radial traversing between inner and outer walls. 
Attached to this traversing apparatus was also a protractor, and with a 
pointer attached to the stem of the probe, angles may be read to within 
l/2°. 


The same constant -cur rent hot-wire set described in reference l4 
was used for turbulence measurements. (A complete description of the 
equipment is given in ref. l6.) Tungsten wires 0.0003 inch in diameter 
and about 0.03 inch long have been employed. The bare tungsten wire 
was copper plated to enable soldering to the probe. The wire was com- 
pensated through a manually controlled resistance-capacitance network 
using the square -wave method. The time constant ranged from 0.7 to 
1.0 millisecond. The mean-square values of the turbulent quantities 
were measured through a vacuum thermocouple unit. Time derivatives of 
the fluctuating quantities were measured with an electronic differentiator. 
All three components of the turbulence intensity, as well as the three 
turbulent shear stresses, were measured by a single inclined wire. The 
technique employed is more fully described in appendix A. In addition, 
an X-wire was used to measure the spectrum of radial fluctuations. All 
spectra measurements were made on a Hewlett-Packard wave analyzer . 

ANALYTICAL CONSIDERATIONS 

With the assumption of axial symmetry, the Reynolds equations in 
the cylindrical coordinate system are 
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V r 


av 9 

dr 


+ v z 


dV 9 + v r ve 
dz r 


i|.( r 2v r ve) - ^(vev,) + V (vev e 



( 2 ) 


V r 





( 3 ) 


where 


2 _ 4 - ^ 4 - ^ 


= 


-.2 r dr n 2 
dr dz 


The equation of continuity for incompressible, axially symmetrical 
flows is 


drV r dV z 

+ = 0 

r dr dz 


(M 


As mentioned previously, the region of flow to be investigated is 
the free- stream region in which V z is practically constant with the 
radius. (This excludes the boundary layers in which V z drops rapidly 
to zero at the inner and outer walls.) In this region of interest, it 
was found that V z increases only slightly with z. The continuity 
equation then gives an estimate of Vr* The experiments indicate that 
V r is smaller than the turbulence intensities but is still of the same 
order of magnitude as the latter. For this reason the terms involving 
V r should, in general, be retained for this type of flow in equations (l) 
to (3). However, the terms involving viscosity v in these equations 


as well as the term Vr 


dV 2 

dr 


in equation (3) can be neglected for the 


region under investigation. (Naturally none of these terms can be neg- 
lected in the boundary layers near the two walls.) 


Although every quantity in equations (l) to ( 3 ) has been measured 
in this investigation, it seems futile to try to balance these equations 
by experimental data. For example, in equation (l) the two terms Vg 2 /r 
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and i — are of a larger order of magnitude when compared with the rest 

P Sr 

of the terms so that any experimental error in these two terms would com- 
pletely overshadow the others. It appears that no useful purpose can be 
served by balancing Reynolds equations except perhaps as a check of 
measurements . 


Among the turbulence quantities to be studied under the influence 
of a gross rotation, the most interesting quantity appears to be v r VQ. 


Idealize the situation a little more and imagine the following case: 

Let the wind tunnel with the inner and outer walls rotating (as a unit) 
at a constant angular speed a) be filled with turbulent flows whose mean 
velocity at the initial instant t = 0 is purely tangential and is such 
that Vq = cor and hence there is no slipping at the walls. The question 
is what will be the V 0 distribution at later instants. The appropriate 
equation describing this idealized case is, corresponding to equation (2) 


but taking into consideration that possibly 



ave / 
at f 


o 


although V r = 


P 


av 9 

at 


i a r 

r2 ar2 


5 2 ( T t + T i) 


( 5 ) 


where for turbulent shear stress 


T t = -pv r v 0 


and for laminar shear stress 


T i = 



It is known that eventually, as the turbulence dies out, the flow becomes 
laminar and, under the condition of rotating walls, the V 0 distribution 

of the fluid must again be that of solid-body rotation, that is, Vq = or. 

It is plausible therefore to assume that, between t = 0 and t = <», 

Vq remains always our. With this assumption, 
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so that either 



T t + T l = 0 


or 


T t + T Z = 


_c_ 

r 2 


The latter possibility must however be ruled out, since it would impose 
a nonzero value for at the walls where t^. vanishes (tj is zero 

everywhere for Vq = our) . 


The statement that v r v 0 is likely to be zero for solid-body rota- 
tion of the mean flow would have been so obvious as to become trivial 
if it were not for the fact that such a conclusion apparently contradicts 
the prediction by the well-known momentum transport theory or turbulent 
shear flows. As pointed out by Taylor (ref. 17) , this theory gives for 
shear flow along circular paths 


r t = pK — — — (rV e ) 
r dr 


(6) 


where pK is the eddy viscosity and tVq is the angular momentum. This 
is analogous to unidirectional shear flows for which = pK — . The 

dy 

basis for equation (6) is, of course, the conservation of angular momen- 
tum of fluid elements in circular flows, which is analogous to the con- 
servation of linear momentum in unidirectional flows. According to equa- 
tion (6), the turbulent shear stress = -pv r vg for solid-body rotation 

of mean flow would be a positive constant and equal to 2pKu> everywhere. 

Now the shortcomings of momentum transport theory are, of course, 
well known. The assumption, for example, that the masses of fluid which 
transport the momentum from one layer to another move about without losing 
their identities and without the influence of viscosity and fluctuating 
pressure gradients is not true. Nevertheless predictions by this theory 
on turbulent shearing stress in unidirectional shear flows are at least 
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qualitatively correct, although quantitatively inexact. The fact that 
in solid-body rotation this theory appears to run into difficulty even 
in a qualitative way is somewhat puzzling. 

Dr. Hsuan Yeh is responsible for the following application of the 
momentum transport theory which makes the theory predict a zero turbulent 
shearing stress for a solid-body rotation type of mean flow. This cor- 
rection is based on taking into consideration the mean pressure gradient 
normal to the flow, a quantity existing prominently in curved mean flows 
but not in straight mean flows. The other simplifying assumptions of the 
transport theory are retained. The detail of such a development is given 
in appendix B. The conclusions can be stated briefly as follows. If at 
any instant a fluid element possesses an excess velocity, either in v r 
or in V0 or both, it will then oscillate around an oval path relative 
to the mean flow. The frequency of oscillation, irrespective of the 
magnitude and direction of the initial disturbance, is always 2cjd. This 
confirms the speculation by Morgan (ref. 12, see "introduction") that 2 od 
might be the resonance frequency of forced oscillations; this value is 
indeed the natural frequency of the system. Furthermore, the amplitude 
of v r is equal to that of vq for any kind of initial disturbance. 

This again seems to confirm Morgan's result for a = 2ca. However, unlike 
Morgan's assumption, the fluctuating velocities v r and vq do not vary 

sinusoidally with time, although the departure from a sinusoidal function 
is not very great. Another interesting fact is that the magnitudes of v r 
and vq at subsequent times do not bear a simple proportion to the magni- 
tude of the initial disturbance. Finally, such a simplified analysis shows 
that when a fluid element moves out with a positive value of v r , for 
example, then there is equal chance of its possessing eith er a positive 
or a negative value of V0, so that the correlation v r v Q is zero. 

The equations governing the six components of Reynolds stress tensor 
(three component intensities and three shear correlations) and the total 
turbulent energy, in cylindrical coordinates, were derived in reference 18. 
These equations, after simplification by the assumption of axial symmetry, 
are reproduced as follows: 



- v r v 0 


2V 0 


r 




v r v 0 2 


r 



dv r 

dz 



v 


( 7 ) 
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V r 


dr 


+ V 7 


dVg2 

dz 


1 d 


dV 0 


+ VrVe ? s; (rVe> + V * T0 5T + 1,0 r = 



lL dq 2 „ dq 2 ) ~ v r — 2 dv r p ^ v z dV 0 

5V' r S?“ + Vz ^j + ve ~ + Tr 5T + v ‘ aT + V * V0 a7 + 


/dVz dV r \ d /V0\ 

V'BT + aT + VrV0r a? T * 


1 d_ 
r dr 


rv T 







/ \ 1 

A v z (£ + ^ 

- V 

dui d Ui ( v 0 2 + v r 2 J 

'v e dv r v r dv e | 

dz ^2 P ) 


d X<j dxj r 2 

\r r d0 r r d0y 


( 10 ) 


rH | C\J 
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(dv r v r 


V 0 ' 


^sr + r/ + v ® r t - 


SV^ v r ^ ^ \ ^Vq 

V9 ST + — ^ rVe) + ^ 5T 


= - — ( Vn ^ + V„ 


2 P\ dr r d9y 


v r 5 


i 3 5“T“ v e^ d Sv e dv r . av e dv r 

2\r2 * ' T + dl We (dT dT + ^ 7* + 


9 . d 



<^e <^r v r v 9 

Sr Sr rC 


(11) 


... ^v z v e ,. ^v z v e \ . l 

r Sr + Z Sz ) 2 


1 d 


V T 


v r v z - — (rVe) + v z v e — + 
r Sr r 


dV 7 


2 Sv 9 


dv„ 


V * V0 ST + Vz " ST + VrVe Sr 


Jjv« 


Sp dp 


Vfl h v, 


2p V 0 dz 


d0 


-i — — r 2 v r v 0 v z + v z ^v 0 


2\ r 2 dr 


_d_ „_2, 

Sz 



_ VI v e av z + !le_ ^Z_ + 
Sz Sz r S0 r S0 


dv e Sv z v 2 v e v z dv r 
+ 


dr dr 2r 2 r r d0 


( 12 ) 
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, Sv r dv z \ 

v r v z l— + — j + v z v 0 


2V 0 ' 


o" aV r — av 2 
v z 2 _ + v r ^ — ± 
oz or 


^ + v. |E - i/l irv r 2 ' 


2p\ dr dz / 2\r dr 


rvr v z - 


VzVe^ . d rr^l ,.l av z ay r . dv z dv r dv z dv r , 

_ + _v rV I . v __ + __ + ___ + 


v r v z v z ^v 0 


2r 2 r r^O 


(15) 


In any one of these equations, the interpretation of the groups of 
terms, from left to right, is as follows: (l) The convective rate of 

change of some turbulent quantity; (2) the production of that quantity 
from the mean velocity (a positive production is associated with a nega- 
tive value for this term); (3) the increase or decrease of the quantity 
due to the fluctuating pressure velocity correlation. This usually serves 
to feed the smaller intensity component at the expense of the larger 
(refs. 19 to 21), as will be more fully discussed later; (4) the increase 
or decrease due to turbulent diffusion; (5) the change due to the action 
of viscosity. This term is a combination of the rate of dissipation of 
turbulent energy to heat and the rate at which viscous forces are doing 
work on a fluid element (see ref. 22). This term is obtained by a trans- 
formation in each equation of the type 


v 



+ v 


■j^ v i 


= -vA^Vj^Vj - 2v 


dvj dvj 
dxi^ dxjj- 


The second derivative of the fluctuating quantity is neglected because 
it is small compared with the other terms except in the region very close 
to the wall. 
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The most direct influence of the rotation would seem to be in the 
production terms, where the peripheral mean velocity appears. The way 
the various turbulent quantities are changed by rotation from the values 
without rotation would depend therefore on this distribution. Since the 
turbulence problem is nonlinear, the mean peripheral velocity of course 
depends in turn on the turbulent quantities. The system of equations 
presented so far is indeterminate. As an equation is derived for each 
of the -unknown fluctuating quantities, new unknowns appear. The theoret- 
ical motivation of deriving such a chain of sets of equations seems to 
be the hope that eventually the new unknowns generated may be evaluated 
by some plausible physical assumption to which the system is not sensitive 
and the process stopped. Meanwhile, the equations have served as a guide 
in interpreting experimental work. In the next sections the determina- 
tion of the mean velocity field and the measurement of some of the various 
turbulent quantities encountered will be presented. 

The following speculation by Dr. Hsuan Yeh on the conditions for simi- 
larity to exist in a turbulence field under gross rotation is of interest. 
For the idealized case in which V r =0, Vq = am, and V z is constant 

everywhere, the equations can be reduced considerably. Since V z now 
serves merely as a carrier, the magnitude of V z cannot influence the 
problem. With respect to a coordinate system moving with a velocity V z , 
equation (7), for example, appears as follows: 


dv r 2 

~dt“ 


2asv r V0 = 


v r dp 
P dr 


1 1 _d_ 

2 r dr 



+ 


v r ve 2 

r 



vr 2 + 2 vr dve 
r 2 r r de 


~ 2 

The last term is proportional to v — where v is a typical turbulence 

A 2 

intensity and A is a typical turbulence microscale. This equation can 
be made dimensionless by dividing each term by dv 2 . ^In the resulting 
equation there appears two dimensionless parameters v/or and v/d\ . 
Physically the first parameter is the ratio of the Coriolis forces due 
to turbulence to the centrifugal force from the gross rotation of the 

v 2 

mean motion. On account of the fact that the dissipation term v — 

A 2 

is proportional to v3/l where Z is a typical length of the energy- 
containing eddies (ref. 23 ), the second parameter can also be written 
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as v/cal and therefore represents the ratio of the Coriolis forces due 
to turbulence to the centrifugal force due to mean rotation of the large 
eddies . 

Since only one rotational speed and one type of turbulence were fully 
investigated in this work, no conclusions can be drawn as to the influence 
of a continuous variation of these two parameters. 


PRELIMINARY MEASUREMENTS 


In order to isolate the effect of rotation on the flow field, it 
was necessary to investigate the flow both with and without rotation. 

As previously described, a prerotation is induced in the flow before it 
enters the rotating screens by a rotating impeller. Due to the twist 
of its blades this impeller when it is not rotating blocks the flow coming 
from the main blower; it must therefore either be modified or removed for 
tests without rotation. 

Initially, the latter arrangement was used. A series of radial 
rods of the same spacing and thickness as the impeller vanes was mounted 
on a dummy hub in an attempt to give similar conditions upstream of the 
system of 20 -mesh screens. These screens were then attached to the hub 
in the same way they had been mounted on the impeller. 

A number of traverses were made at various axial positions both with 
and without rotation, using the dummy hub in the latter case. The speed 
of rotation of the tunnel was 1,560 rpm (corresponding to maximum rated 
speed of the driving motor), giving a peripheral velocity of 57*7 fps 
and 72.8 fps to the inner and outer walls, respectively. A throughflow 
of 19.5 cu ft /sec was used both with and without rotation, corresponding 
to an average throughflow velocity of 42-5 fps. This value was arbitrarily 
chosen to give a ratio of peripheral to axial velocity of nearly unity 
near the inner wall. The Reynolds number based on this velocity and half 
the distance between inner and outer wall is 26,700. 

The results of these measurements are shown in figures 6, "], and 8, 
where axial velocities without rotation, axial velocities with rotation, 
and peripheral velocities, respectively, are shown at various axial sta- 
tions. The axial distance z is measured from the farthest downstream 
20-mesh screen. (it should be noted that the 4 -mesh turbulence -producing 
grids mentioned in the section "Wind Tunnel" were not yet installed.) 

The axial velocities are made dimensionless with the average throughflow 
velocity, the peripheral with the peripheral velocity of the outer wall. 

In figure 8, the striking lack of change of the peripheral velocity 
with distance along the axis should be noted. Although the peripheral 
velocity is not strictly one of solid-body rotation, it nevertheless has 
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changed very little when compared with the changes that have occurred in 
the axial velocity profile in the same distance downstream. 


Turning to the axial velocity profiles as a function of distance 
downstream both with and without rotation as shown in figures 6 anri 
it is interesting to compare the thickness of the boundary layers on 
both inside and outside walls for these profiles. The boundary- layer 
thickness is here defined as that distance from the wall at which the 
velocity distribution given by the power law which best fits the observed 
data is equal to the free-stream value. Results are given in figure 9 . 

It is seen that without rotation the boundary layer on the outer cylinder 
is thicker than that on the inner. For fully developed laminar flow in 
an annular space the maximum velocity occurs closer to the inner wall 
than to the outer wall (see ref. 24). To satisfy this condition the 
boundary layer on the outer wall in a laminar flow that is not fully 
developed will have to grow faster than that on the inner wall. It is 
interesting to note that the same trend still holds for turbulent flow. 


With rotation, the very interesting effect is seen that the boundary 
layer on the inner wall grows faster than it did before, while the oppo- 
site is true for the boundary layer on the outer wall. A detailed analysis 
of boundary- layer growth must include measurements very close to the wall. 
This is beyond the scope of the present report, since the emphasis is on 
the region between the boundary layers . 

A profile of turbulence intensity v z /V z at the farthest downstream 
station without rotation is shown in figure 10. As would be expected, 
at a given distance from the wall a larger turbulence level exists near 
the outer wall than near the inner, since the boundary layer near the 
outer wall is thicker. 


From the 
of magnitudes 


velocity profiles just presented, one may estimate the order 
of the following terms in Reynolds equations (l) to ( 3 ): 


v r 


dV r 

dr 
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ESTABLISHMENT OF SIMILAR UPSTREAM CONDITIONS 


The measurements and results described in the previous section were 
obtained, in the case of nonrotating flow, with the impeller removed from 
the system and replaced by the dummy hub and rods described. This was 
intended to give the same conditions upstream of the screens both with 
and without rotation. The effectiveness of this device cannot rigorously 
be established without a direct comparison to nonrotating flow produced 
with the impeller left in place . 

It was found that if the impeller is rotated at a low speed, a flow 
with no peripheral velocity is obtained behind it. This is so because 
the trailing edge of the impeller vanes is not axial but rather at an 
angle of about 20° to the axis of rotation. The flow is deflected toward 
the axis by the 20-mesh screens behind the impeller but not enough to 
make it social. A slight rotation (about l60 rpm) then eliminates the 
remaining peripheral component. At this low speed of rotation, however, 
the inlet vane angles are not correct. An upstream extension was there- 
fore attached to each moving blade , tangent to it at its leading edge 
and designed to meet the inlet flow at the proper angle. To reduce 
secondary flows between impeller vanes, a series of four flat concentric 
rings was attached to the vanes by cutting slots chordwise into the 
trailing edge of the vanes and fitting the rings into these slots. These 
rings, l/2 inch wide and coaxial with the axis of rotation, extend 
l/4 inch upstream and l/4 inch downstream of the vane trailing edge. The 
20-mesh screens were reattached. A mean velocity profile was then 

obtained for slow rotation at — - — = 1.1, corresponding to the farthest 

b - a 

upstream profile shown in figure 6. At this station, comparison of the 
profiles with slowly rotating impeller and with dummy hub and rods showed 
them to be almost identical. 

Preliminary turbulence intensity measurements in the mean flow direc- 
tion showed the turbulence level also to be about the same with the two 
different methods of producing purely axial flow. Thus the data taken 
with the dummy hub with its radial rods shown in figure 6 are believed 
to offer a good comparison with the flow in rotation. 

With slow rotation the axial turbulence intensity component at 

5 — = 3.2, a traversing station fairly close to the screens, varied 

b - a 

in peripheral direction from 1.8 to 2.2 percent of the mean velocity with 
the impeller in place. These variations were traced to wakes from the 
impeller vanes passing through the four 20-mesh screens. At this dis- 
tance from the screen, which corresponds to about 150-mesh lengths, the 
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turbulence created by these screens has therefore decayed to such a small 
value that nonuniformities passing through the screen become significant. 
In order to avoid these axially asymmetric blade wakes one could go 
farther downstream until they have died out, but one then approaches the 
flow region where the axial velocity is no longer uniform. It was decided 
to install additional large-mesh screening some distance downstream of 
the 20-mesh screens. A large-mesh screen will create turbulence of a 
larger scale which allows more time to decay. At a given station, the 
turbulence level will then be higher, and if the peripheral variations 
in the turbulence upstream of the screen remain the same, such variations 
downstream of the screen will be proportionately lower. The position of 
the screen will be a compromise: If it is too fair downstream the boundary 

layers are too thick, too far upstream the wakes become stronger. 

Various arrangements of 4-mesh, 0.047- inch-diameter wire (Solidity = 
0.340) screens were tried. The best arrangement turned out to be the 
following: one 4-mesh screen l60 mesh lengths (referred to the fine 

screen) behind the last of the 20-mesh screens followed by another 4-mesh 
screen 23 meshes (referred to the coarse screen) farther downstream. The 
traversing station closest to this screen is 17*5 mesh lengths away. 

(This station will henceforth be identified by = 17-5 where I is 

the axial distance measured from the last 4-mesh screen (in contrast to 
z which is measured from the farthest downstream 20-mesh screen) . All 
further radial traverses were made at this station. In terms of the 
distance between inside and outside walls, this station corresponds to 

—2 — = 7.1. Peripheral variations at this station were only ±2*5 percent 
b - a 

of the turbulence level. Measurements taken at this station are presented 
in the following section. 


EXPERIMENTAL RESULTS 


Mean Velocities 


Axial velocity profiles with and without rotation, 17*5 mesh lengths 
downstream of the last 4-mesh screen, are shown in figure 11. The through- 
flow velocity is the same as that for the preliminary experiments, but 
with rotation the angular velocity has been increased from 1,560 to 
2,000 rpm to maximize the effects of the rotation within the limits of 
the equipment. It can be seen that the presence of the two additional 
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screens has made the axial velocity more uniform compared with the pre- 
vious profiles shown in figures 6 and 7* As before, the rotation seems 
to have decreased the boundary layer on the outside wall, which without 
rotation is thicker on the outer wall than on the inner. The peripheral 
velocity is shown in figure 12, where for comparison the previous profile 
is shown as a solid line. Due to the higher speed of rotation and the 
two extra screens, the peripheral velocity distribution is much closer 
to solid-body rotation than it was before. Figure 12 shows that there 
is still a small relative peripheral velocity between fluid and screen. 


Turbulence Intensity and Shearing Stresses 


The method used in measuring the various fluctuating quantities is 
described fully in appendix A. The essential point is that with this 
method the three components of the intensity as well as the three shear 
stresses are measured at one time. Results of these measurements follow: 


Figure 13 shows the radial, tangential, and axial turbulence inten- 
sity, without rotation, made dimensionless with the local mean velocity. 
As would be expected for a grid-produced turbulence with no mean velocity 
gradients, the intensities are fairly uniform in the region excluding the 
boundary layers, with Vr/V z and vq/V z about the same, both less than 
^z/Vz* Again, because of the thicker outer boundary layer, the intensity 
near the outer wall is larger than that near the inner wall, at the same 
distance from the walls. The three shear correlation coefficients for 
purely axial flow are shown in figure 14. Since there are no mean veloc- 
ity gradients in this case except for those near the wall, one would 
expect no shear stresses other than those near* the wall. Figure 14 shows 
that this is indeed the case. 


For experiments with rotation, two sets of results are presented. 

One set is in a s,n,r coordinate system, with s in the mean velocity 
direction, n perpendicular to it in the stream surface, and r radial. 
The other set of results is in the usual cylindrical coordinate system. 
The conversion from one to the other is given in appendix A. 

With rotation, figure 15 shows the three components of the turbu- 
lence intensity in s,n,r coordinates (made dimensionless with the local 
mean velocity) . It is to be noted that, except near the outer wall, 
v s /V s and v n /V s are nearly the same. The radial fluctuation intensity 
is much less. Shear correlation coefficients in the same coordinates are 
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given in figure l6; although there is considerable scatter of the data 
it appears that the coefficients are no longer zero. 


The same results in r,0,z coordinates are shown in the next two 
illustrations. The intensities are given in figure 17 and were made 
dimensionless with the axial component of the mean velocity. Since the 
axial velocity is nearly the same both with and without rotation, the 
results as plotted can be directly compared with the intensities pre- 
viously given for purely axial flow. The rotation has not appreciably 
changed the axial fluctuations, but the radial fluctuations have been 
markedly decreased and the peripheral, increased. The shear correlation 
coefficients are shown in figure 18. The coefficients R rz and R z q 
still appear to be zero, the same as they are with no rotation. How- 
ever, the rotation is seen to have produced a slightly positive value 
for R r 0 over the entire flow area. Since in the peripheral Reynolds 

equation, equation (2), the term — (r^VpVg) appears, the variation of 

dr 


in figure 19, with V Q again the 

average throughflow velocity. 


v r v 0 


is shown as 


v r v 0 


b - a 


A more direct comparison of the turbulence intensities with and 
without rotation is given in figures 20, 21, and 22. 


Integral Scales 

The integral scale in the direction of the mean velocity is defined 
as 



dS 


where R s is the correlation between simultaneous measurements of v s 

at two points separated by the distance s, measured in the mean flow 
direction. From the ratio of the compensated to uncompensated signal 
from a hot wire placed normal to the flow, Ls may be determined approxi 
mately (see, for instance, ref. l) . Results for Ls/M, with M the mesh 
size of the screen, are shown in figure 23 both with and without rotation 
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As a comparison, reference 26 shows = 0.22 for the same ratio of 
distance to mesh size behind a single 1-mesh screen at 32.8 fps. 


Microscales 

The microscales of turbulence are here defined in the following way: 



The mean- square derivatives appearing above can be related, in 
homogeneous turbulence, to the various correlations of the particular 
velocity component at a point with the same component at another point 
removed from it in the proper direction. (See, for instance, ref. 1.) 

In isotropic turbulence all these microscales are the same, but this is 
not true otherwise. The measurement of two-point correlations was not 
made in the present work. However, by Taylor's hypothesis, ds may be 
replaced by -V s dt, where ds is in the mean flow direction and V s 

is the mean velocity. Hence, a direct measurement of the time derivative 
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may be used to calculate microscales, but only those involving derivatives 
in the mean flow direction. 

For flow with no rotation, only one radial station was measured at 

= 0.5, halfway between the two walls. Since the turbulence level 

b - a 

was found, to be fairly uniform in the region away from the boundary layers, 
it was not considered necessary to take a complete radial traverse. It 
was found that, in this region, 


A v z ,z 


M 


0.17 


^v fl ,z ^v , z 

= — £L- = 0.19 

M M 


For the rotating case, since only the derivatives involving the 
streamwise direction s can be measured by time differentiation, scales 
in the streamwise direction are given in figure 24. To express these 
scales in the more meaningful cylindrical coordinate system, it becomes 
necessary to make certain assumptions based on isotropic relationships. 

The details of such transformation of coordinates are given in appendix A. 

^v z ,z ^v 9 ,z ^v r ,z 

With these assumptions, , , and were calculated. 

MM M 

They are shown in figure 25* The values without rotation in the region 
excluding the boundary layer are indicated in the same figure by dashed 
lines . 


Turbulence Decay 

Since the flow is not fully developed, it is important to know how 
the turbulence is changing with distance along the axis. All of the 
measurements reported so far in this section have been made at a fixed 

distance from the farthest downstream screen, at - 7 = 17*5* With suit- 

M 

able modifications in the traversing apparatus to allow movement in the 
axial direction, axial traverses were made over a distance of about 
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10 mesh lengths at the midstream position without rotation and at several 
radii with rotation. 


Without rotation, measurements of v z /V z agreed very well with 
some screen decay data by Corrsin (ref. 26) 


VQ 

Since — 

V z 


— at 

V z 


t 


- = 17 - 

M 


was found also to be in agreement with this reference, the variation of 
these latter quantities in the axial direction was taken from there. 
Figure 26 shows the inverse of the turbulence levels plotted against 
distance from the screen, again made dimensionless with the mesh size. 


As is described in appendix A, the method of measurement of the tur- 
bulence employed herein makes use of the signals from a single hot wire 
oriented in various positions with respect to the mean velocity. In each 
position, it is sensitive to one or more of the turbulent intensities and 
shear stresses; by an algebraic manipulation of the signals for the vari- 
ous positions the desired quantities are obtained. It was found that with 
rotation of the flow the wire signals as a function of axial distance 
could be collapsed into a single curve, regardless of radius or wire posi- 
tion. This implies that for the range of axial distance considered, the 
ratio of any one intensity component to another remains constant, within 
the accuracy of measurement. Making use of this fact, decay curves with 
rotation are shown for the three intensities at the midstream position 
in figure 27. For any other radial position, the curves are shifted in 
such a way as to go through the previously measured values at that posi- 
tion, at ^ = 17.5 • It can be seen that the effect of rotation has been 
M 

to lessen the decay, for all components, with the strongest effect on the 
radial component. 


Spectrum- 1 

The spectra of v r 2, vq 2, and v z 2 with rotation are shown in 


figure 28, all taken at one point 


= 0.2, 1 = 17.5 

ib - a M > 


The abscissa 


of these curves is the one-dimensional wave number Kj_, oriented along 
the direction of the mean velocity. In figure 29 the v r 2 spectra at 


two radial positions, 


= 0.2 and 0.8, are compared. 


1 


Spectral measurements were made by Dr. Hsuan Yeh. 
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DISCUSSION 


General Considerations 

The left-hand sides of equations ( 7 ) to (13) for turbulence inten- 
sity and shear stress components become, respectively, when simplified 
by Vq = oar in the present problem, 



v r VQ2oa 


1 dv e 

- V z — — + v r vg2oa 

2 dz 




d Y e v z 

dz 


+ 


v r v z a) 



^ v r v z 

6z 




The statement on orders of magnitude (given in the section 'Prelim- 
inary Measurements") for the various mean velocity terms has been used 
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to eliminate those terms that would otherwise appear. The terms on the 
right-hand sides of equations (7) to (13) remain the same and are not 
repeated here. 

In the simplified left-hand-side expressions just given, the second 
term stands for the production of the quantity whose rate of change is 
given by the first term. In the present case, this first term is found 
to be negative for the first three equations, corresponding to a decay 
of turbulence . A negative sign for the second term implies a positive 
production, and vice versa. For example, in the first expression, if 

v„vq < 0, dv 2 /dz would be decreased by the shear stress acting on the 

1 — / 

mean velocity. If v r Vg *-*, 9 v r Ibz would be increased. A pos itiv e 

production of course does not necessarily mean a larger value of v r 2 
compared with no production since the pressure and triple velocity cor- 
relation terms also enter into the energy balance. 

The influence of rotation is directly apparent in the production 
terms; with o> = 0 all these terms vanish. However, all the terms on 
the right-hand side of equations (7) to (13) will probably also be dif- 
ferent with and without rotation. 

Flow without rotation .- In considering flow without rotation, all 
production terms vanish, and the turbulence decay is given then by the 
combined influence of the viscous terms, pressure correlations, and triple 
velocity correlations. At any point the values of the turbulence inten- 
sities depend on this decay and on the manner of production at the screen. 
Immediately downstream of the screen, velocity gradients exist in the 
flow due to the wire wakes. There are both dV 2 /r bd and bV z /br terms, 

alternating in sign on either side of a wire wake, (it would be more 
natural in the nonrotating case to use Cartesian coordinates, with x,y 
instead of r,6. The cylindrical coordinates are retained, however, for 
comparison with the rotating case . ) 

Very close to the screen, the left-hand sides of equations (7) to 
(9) and (ll) to (13) become for turbulence intensities: 


1 

2 



dz 




dz 
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v 8 v z 


av z 

r Se 


+ 


V V 

r z 


8v z 

dr 


and for turbulent shear stresses: 


1 dv r v 9 

2 Z dz 


1 y d v 6 v z [ v 9 2 8v z | v r v 9 8 v 2 

2 Z dz 2 r ^9 2 


l v dv r v z , v r^ dV z , v r v 9 dV z 
2 Z dz 2 Sr 2 r d0 


Assume that, in the shear equations only, the production term alone 
is important in the determination of the value of the quantity being pro- 
duced. By this assumption, since there is no production of v r vg, this 

quantity vanishes. In the last two shear-stress equations, it follows from 


the same assumption that when 


av z 

r Se 


> 0 


then 


v e v z 


< o 


and when 


< 0 then v Q v >0. It also follows that when > 0 then 


r d0 


'9 v z 


Sr 


v r v z < 0 and when 


Sv z 

Sr 


< 0 then v r v z > 0 


Thus the quantities v 9 v z SV z /r S0 and v r v z SVz/Sr are negative 

throughout the entire region immediately behind the screen. Co nsid ering 
now the intensity equations, there is no direct production of Vg2 and 

v r 2 but always a direction production for v 2 ^ • It is known that some 
distance behind the screen the turbulence will become nearly isotropic. 
Since there is direct production from the mean velocity only of v z 2, 
the other two components must receive their energy from the pressure 
velocity correlation (see, e.g., refs. 19 to 21) . Thus the assump- 
tion made above for the shear equations that the production term alone 
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determines the value of the quantity being produced cannot be valid for 
the intensity equations. 

p 

The physical picture is then a direct production of v z , with about 
equal feeding of Vq 2 and v r ^ by the pressure velocity correlations. 

Some distance behind the screen when the wakes have died out there is no 
longer any direct production since the velocity gradients in V z have 

vanished; the terms v r v z and VqV z therefore vanish, while v z 2 , 

Vq 2 , and v r 2 decay, with v z 2 larger than Vq 2 « v p 2 but decaying 
faster due to energy loss by pressure fluctuations which feed Vq 2 
and v r 2 . 

Flow with rotation .- A similar reasoning is now applied to the case 
when rotation is imposed on the flow field. Far enough downstream so 
that the wire wakes have disappeared, the only production terms are those 
given in the section "General Considerations.” 

Within the limits of the previous assumption it is seen that in the 
last two equations the production of VgVg necessitates a nonzero value 

of v v , while the production of v v necessitates a nonzero value of 

r z r z 

VgVg. Thus v r v z is necessary to produce VgV z which is necessary to 
produce v r v z and if v r v z is zero it must remain so. The same is 
true for V 0 V z . Furthermore, the system is such that if in the vgv z 
equation, v r v z > 0 , the production of vgv z is negative, and vgv z < 0 . 
Then in the v r v z equation the production of v r v z is negative so that 
v r v z decreases. Arguing from the production terms alone, the last two 
equations are such that neither vgv z nor v r v z could be produced 
directly from the rotation. 

The equation for v r vg, however, involves the turbulent intensities 
in the production term. It is seen that if v r 2 > V q 2 then v r vg < 0 
and if v r 2 < Vg 2 then v r vg > 0. 

Thus if v r 2 > vq 2 , vg 2 is produced from the rotation by the 
turbulence- intensity equation, while v r 2 is diminished. Similarly, 
if vq 2 > v r 2 , vq 2 is diminished while v r 2 is produced. When 
v r 2 = vq 2 , there will be no production of v r vg and thus no effect of 
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rotation on the production of either ■ v r 2 or vq 2. Hence the production 

terms act in such a way as to keep v r 2 « v 0 2, even with rotation. The 

actual values of v r 2 and v 0 2, of course, will depend on the more 

indirect effects of rotation on the remaining terms in the first three 
equations . 


Close to the screen, the production of turbulence is not affected 
by the rotation if there is no relative motion between screen and fluid, 
that is if both axe in solid-body rotation. If there is relative motion, 
however, similar reasoning as before shows a positive production of 
dV 0 

v q 2 due to the term v z ve — -, which is found to be always of one sign. 

dz 


'9 


The physical picture, then, for solid-body rotation is a direct pro- 
duction at the screen of v z 2 , with subsequent direct production or dimin- 
ishing of v r 2 and ve2, depending on the values of v r 2 and vq2, in 

such a way as to keep them about equal. If there is relative motion 

between fluid and screen there will be an additional production of v 0 2 

at the screen. The turbulence then decays under the influence of viscous 
dissipation and energy transfer between components by the pressure velocity 
correlations; both quantities are presumably also influenced by the 
rotation. 


Energy Balance 

In considering the energy balance it will be recalled that, at the 
one axial station considered, the experiments showed the following 
results over most of the flow field, excluding the boundary layers: 

Without rotation, 


v z 2 > v r 2 « v 0 2 


v r v z = v 9 v z = v r v 9 = 0 
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with rotation, 


v z 2 « v 0 2 > v r 2 
= v e v z = °> v r v 9 > 0 


The result that v r vg >0 is entirely consistent with the finding 

of v 0 2 > v r 2 by the preceding considerations. Since measurements of 

the microscales and turbulence decay were made in addition to the meas- 
urements of the quantities above, it is possible to substitute these 
measured quantities into the three turbulence- intensity equations (7) 
to (9). 

For the purposes of calculation, the turbulence -intensity equations 
were made dimensionless by multiplying each term' of these equations by 

m/v o ^. In calculating the dissipation term, the following assumptions, 
which are strictly valid only for isotropic turbulence, were used: 



An assumption of this type is necessary because only derivatives in the 
flow direction were measured. Furthermore , since neither the triple veloc- 
ity nor the pressure velocity correlation terms are measured, they will be 
considered together. Substituting the measured values into the energy equa- 

£ 

tions for the case of no rotation, it was found that, at -- = 17*5 and 

M 


r - a 
b - a 


0.5, (without rotation, the energy balance is made only at the 
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midstream position, since conditions were previously found to be radially 
■uniform except near the wall) : 



Decay 

term 

Production 

term 

Pressure velocity and triple 
velocity correlations 

Viscous 

term 

v r 

5.6 X 10 “ 5 

0 

-0.3 x 10-5 

3-9 X 10“ 5 

v e 

3.6 x 10"5 

0 

- 0.3 X 10 “ 5 

3.9 x 10“5 

v z 

6.4 x 10~ 5 

0 

0.6 X io -5 

5.8 X 10-5 

q 

13.6 x 10 “5 

0 

0 

13.6 x 10"5 


The exact balancing of the total energy equation is not a coin- 
cidence. Measured values of the derivative terms consistently gave 
A too large and a dissipation which was too small to balance this 
equation. This discrepancy, amounting to a factor of about 2.0, could 
be due to the assumption of isotropy used in their calculation. For 
lack of a better procedure, all measured values of A were adjusted by 
the same amount in such a way as to balance the total energy equation 
for the case of no rotation. 

A result consistent with the previous considerations is obtained. 

The unknown pressure-velocity-correlation and triple-velocity-correlation 
term is obtained by difference; when positive it represents a loss of 

energy and when negative, a gain. Thus v r 2 and vq2 are seen to be 

gaining energy at the expense of v z “^ due to the influence of this 
effect . 

The same procedure when applied at several radial positions in the 

rotating flow at the same station, i = 17.5, gives the results shown 

M 

in figures 30 to 33 for the three intensity components and the total 
energy. In each case, the relative gain and loss due to the various 
terms necessary to make up the measured decay are shown. The viscous 
term, combining dissipation and viscous work, is denoted by DVW. The 
pressure velocity and triple velocity term, considered together and 
denoted by PTC, is again obtained by difference. 
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The balance on total energy shows over the center region a larger 
loss due to viscous action than the decay, which implies a gain due to 
the term obtained by difference. Near the walls the opposite is true. 
This implies a flow of total energy from both walls to the center. The 

same considerations hold for the energy balance on v z 2 . Similarly, for 

both v 0 2 and v r 2 , a comparison of the decay term with the viscous term 

shows the former to be less in the center, implying an inward flow of 
energy. The picture is made more complicated by the production term 

2v r v 0 cD, which for v r VQ > 0 gives a gain by direct production for v r 2 

and a loss for Vq 2 . This must be made up by the pressure and triple 

correlation term, which therefore gives a gain for v 0 2 and a loss for 

v r 2 . It will be recalled that this term, without rotation, acted as a 

P P 

gain for both v r and v 0 • 

It is interesting to compare the coefficients A both with and 
without rotation. These coefficients are defined in the following way: 


A V Z ,Z _ 


M 

& 

>* 

tsi 

1 

A r 

M 
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For a given value of v, an increase in A implies an increase in A 
and therefore a decrease in viscous dissipation and work. Figure 
shows the distributions of A z , k r , and Aq calculated from the 
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measured intensities and time derivatives. It can be seen that A z and 
Aq are about the same, with A r considerably less over the entire range. 

Thus, the relative loss of energy by viscous action is larger for v r 2 
than for the other two components. 

A tentative physical picture for the effect of rotation is then as 
follows: The rotation increases the loss through viscous action for the 

radial fluctuations relative to the peripheral and axial fluctuations. 

A shear stress v r vg is produced by the difference of vq 2 and v r 2 

in such a way as to restore their equality by direct production from the 
mean peripheral velocity. Because of this, the pressure velocity and 

triple velocity correlations, which without rotation fed v r ^ and vq 2 
from v z 2, now feed vq 2 at the expense of v r 2. 

It is interesting to note that solid-body rotation causes an increased 
loss of energy by viscous action for the radial velocity fluctuations com- 
pared with the radial and peripheral fluctuations. It will be recalled 
that in the Introduction reference was made to investigations of small 
motions relative to solid-body rotation. There it is found that the 
radial velocity component is suppressed by the rotation, but this is an 
effect associated with potential (Coriolis and centrifugal) forces and 
not with viscosity. 

It should be pointed out that the considerations above are beclouded 
by the following facts: 

(1) The turbulence -producing screen has a peripheral velocity rela- 
tive to the fluid since the fluid is not quite in solid-body rotation 

at the screen. This implies a larger vq 2 production than would be due 
to purely solid-body rotation, as previously mentioned. 

(2) Due to the presence of the walls and axial velocity boundary 
layers, there is a flow of turbulent energy inward from the walls. 

(3) Because only time derivatives were measured, the viscous terms 
of the energy equations are calculated with the assumption of isotropic 
relations between the derivatives of any one velocity component with 
respect to different directions. This assumption is usually verified 
for sufficiently high values (of the order 100) of R^, defined as 
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For the different intensity components, the values of obtained 

in this investigation are shown in figure 35, with and without rotation. 
It is seen that they are comparatively low, casting some doubt on the 
validity of the assumption of isotropic relations between viscous terms. 


Production of Turbulent Shear Stresses 

The assumption that in the turbulent-shear-stress equations only 
the production term is important in the qualitative determination of 
these stresses will now be further examined. 

It will be recalled that in any one of these equations it was 
assumed that a positive production would give a positive shear stress, 
negative production, a negative shear stress, and no production would 
lead to vanishing of the shear stress. For unidirectional shear flows 
this assumption is certainly true. 

Take the case of fully developed flow through a two-dimensional 
straight channel, where x is the direction of mean flow and y is 
the distance from the wall. For this case, the production terms of the 
shear-stress equations are as follows: 


For v x v z : 


v y v z dV x 

2 dy 


For 


V v • 
x v y* 


v 2 


dV. 


2 dy 


x 


For v v_ : 0 

y L 

Again a negative sign of these terms implies a positive production. Now 

since there is no production of VyV z , v y v z vanishes. There is then no 

production of v Y v_ , which also vanishes. Since v 2 is always > 0, 

there will be a production of negative v x v y when dV x ^c)y is positive, 

and vice versa. Hence v v is negative in the region where dV /dy is 

x y x / 

positive, and vice versa. At the center of the channel ^ = 0, there is 

oy 

no production of v x v y ant ^ v x v y will be zero there. Thus the assumption 
is in complete agreement with the experimental fact and also with the well- 
known momentum transfer theory. 

The case of fully developed flow in a two-dimensional curved channel 
is considered next. An investigation of this flow configuration is 
reported in reference 14. Using cylindrical coordinates, let the 0 
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direction be in the direction of mean flow and r, the radial direction 
(i.e., the direction across the channel). The production terms of the 
shear-stress equations are as follows: 


For 


W 


v r v z / 5v9 + v©\ 

2 \dr r / 


For v v n : 
r 9 




ve 

r 


For 


v v : 
r z 


- v e v z 


Vq 

r 


For a very large radius of curvature, the terms involving Ve/r 

become small, and in the limit of zero curvature the production terms 
are the same as in the straight channel except for the change of nota- 
tion from y to r and from x to 9. The influence of curvature is 
thus represented by the terms involving Vg/r. For a fully developed 

vanishes (this can be shown from the Reynolds equations, 

it also is 


v r v z 


flow, 

see ref. 14 ) ; then, since there is no production of 


zero. Close to the wall, 


v 0 v z> 


dV 


9 


dr 


'9 


» — , and the first term alone produces 


a negative shear v r v g on "the inner wall and a positive term v r v Q on 
the outer, by the assumptions made herein. The influence of curvature 


on v r ve is seen to depend on the term 



If 


'9 


> 


2 v / r o 2 

over the entire flow, this term will cause an additional production of 
everywhere. Thus the v r VQ shear-stress distribution 

over the entire channel 
distributions 


positive v r vg everywhere. Thus the v r vg 
would shift in the direction of positive v r V0 
under the influence of a mild curvature. Measured 


V r V0 


obtained in such a curved channel are reproduced in fig ure 36; the results 
were taken from references 14 and 16. For this flow, vq^ does exceed 


V r ‘ 




and a shift of v r V0 in the proper direction is in fact observed. 


A physical interpretation of the influence of curvature was first 
suggested by Rayleigh (ref. 27) in the following way: It may be assumed 

that for flow with curvature a displaced fluid element conserves its 
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angular momentum. Then a positive radial displacement of a fluid element 
is associated with a velocity defect if the mean peripheral velocity 
decreases more slowly with radius than l/r. The radial motion of the 
element is retarded since the radial pressure gradient existing in the 
mean flow exceeds the centrifugal force acting on the fluid element. On 
the other hand, if the mean peripheral velocity decreases with radius 
faster than l/r, any radial displacement would be increased by the action 
of the radial pressure gradient. Thus radial motion is suppressed on 
the inner wall of the curved channel and increased on the outer, 
resulting in an increase of v r v 0 on the outside and a decrease on the 
inside. 

As is shown in the "introduction'* and appendix B, such simple rea- 
soning requires an extension in the case of solid-body rotation. The 
modified theory still assumes the conservation of angular momentum (which 
is, of course, not exactly true) but traces through the complete history 
of a displaced particle. The theory then predicts that v r v 0 = 0 for 

solid- body rotation. This still does not explain the experimental fact 
that v r Vg is slightly on the positive side. A simple physical inter- 
pretation of the positive value of v r v 0 found in solid-body rotation 
might be as follows: For a fluid element a positive value of v 0 cor- 

responds to a small motion relative to the rotation in the direction of 
rotation. This would be associated with a Coriolis force .acting radially 
outward, causing a v r > 0. Thus v r v 0 > 0. Similarly, a small negative 

peripheral velocity relative to the rotation would give rise to a nega- 
tive value of v r , again giving v r v 0 > 0. Thus a peripheral fluctua- 
tion causes a positive value of v r v 0 . Similar reasoning shows that a 

radial fluctuation c auses , again due to the action of Coriolis forces, 
a negative value of v r v 0 * If the peripheral fluctuations are larger, 

as they are here, the former effect predominates to give a resultant 
positive value of v r v 0 . 


CONCLUDING REMARKS 


A peripheral velocity distribution like that of a rotating solid 
body imposed on the turbulence produced by screens is found to undergo 
very little change with distance along the axis of rotation. Experimental 
measurements of the turbulent shear 1 correlation between radial and periph- 
eral fluctuations 17*5 mesh lengths from the farthest downstream screen 
give a small positive value for this quantity over the major portion of 
the flow field. The other shear correlations are found to vanish. It 
is believed that the nonvanishing of this shear correlation is due to 
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the fact that the flow is not fully developed, although the exact mech- 
anism is not clear. 

Consideration of the qualitative nature of the turbulent shear 
stresses, based on the production term in the equations for these stres- 
ses, leads to results in agreement with the experimental findings. Such 
considerations also lead to results in agreement with experiment for the 
influence of curvature on a fully developed channel flow. 

The influence of rotation on the intensity of the fluctuations 
I 7.5 mesh lengths from the screen is as follows: The radial fluctuation 

intensity is suppressed, the peripheral intensity is increased, and the 
axial intensity is not largely affected by the rotation. The turbulent 
energy balance for these intensities indicates a larger loss of energy 
by dissipation and viscous work for the radial intensity than for the 
other two, with a gain by production from the mean velocity for the 
radial fluctuations, and a loss by production for the peripheral fluctua- 
tions. The action of the pressure velocity and triple velocity correla- 
tions, which without rotation serve to increase the radial and peripheral 
fluctuations at the expense of the axial, is, with rotation, to increase 
the peripheral fluctuations at the expense of the radial. These results 
are obtained under the assumption of certain isotropic relations for the 
viscous terms of the energy equations. 

The influence of rotation on the various mean-square derivatives of 
the turbulent fluctuations which make up this viscous term should be 
investigated further; this would serve to check the validity of the assump- 
tion just mentioned and to aid in the interpretation of the result that 
radial fluctuations experience a relatively larger energy loss through 
viscosity than do the peripheral and axial f luctuations . 

It would also be extremely interesting to study the flow considered 
herein after it has become fully developed. The measured lack of change 
of the peripheral velocity profile with axial distance indicates that the 
length of tunnel necessary to achieve this would probably not be larger 
than that necessary to achieve fully developed flow without rotation. 

The complications introduced by changing conditions along the axis would 
then be eliminated to make a more detailed analysis possible. 


The Johns Hopkins University, 

Baltimore, Md., October 8 , 1958* 
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APPENDIX A 


METHOD OF MEASUREMENTS 


Response of a Hot Wire 

The response of an electrically heated fine wire (the hot wire) in 
a steady flow of air stream is found experimentally to obey the following 
approximate relation when all quantities have reached an equilibrium state 
with time: 


where 



(Al) 


I wire heating current 

R Q wire resistance with no heating 

R wire resistance (with heating) at operating condition 

Vp velocity component perpendicular to axis of wire 

A,B constants for a given wire and surrounding medium, to he deter- 

mined experimentally 

•p 

In practice, — ^ is fixed by the experimenter. In the present 

p O 

investigation is chosen to be 3*0 at all times, corresponding 

R - R 

R _ R ^ -q 

to 2., the overheat ratio, of 0.5. For a fixed value of - — K , 

R 0 H - R 0 

it is convenient to rewrite equation (Al) in the following form: 


I 


2 



(A2) 


2 

The two constants I Q and S are determined experimentally by meas- 
uring over a suitable range of Vp in a calibrating air stream. The 

term I Q ^ has the significance of being the extrapolated value of 


at zero Vp. It is clear that the constants I Q and S are related to 
A and B by I 0 2 = B^ R R °^ and S = a( R -~ — -) . 
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If velocity fluctuations exist in the airstream, the wire tempera- 
ture and therefore its resistance will also fluctuate with time. In 
constant-current operation, I = Constant, the resistance fluctuation dR 
is measured as a voltage fluctuation dE. Differentiating both sides of 
the steady-state equation (Al), in which R and V p are the variables, 

gives 


dVp = 2 I 2 dE 

w -p -r 2 -r 2 2 

V P JL - 1 1 ~ l o 

R o 


In practice, since I does vary slightly, one may have to take its 
variations into account also, resulting in the extra factor called the 
"degeneracy" or 


dVp 



Oi 

M 

CVJ 

1 + 2\ 

fa - Ro) 

il_ 

V R o l 

e B 

1 _ i x 2 _ i 2 

R ± O 

1 

_ JL 

e b 



(A3) 


where Eg is the supply voltage to the Wheatstone bridge in which the 
wire is placed. These equations are valid if the fluctuation occurs 
very slowly (quasi-steady condition) or if the effect of the finite heat 
capacity of an actual wire has been electronically compensated. The pro- 
cedure so far is quite standard; see, for instance, Kovasznay (ref. 28). 
By taking the root mean square of equation (A3) , the root -mean-square 
value of the velocity fluctuation is obtained in terms of the root -mean- 
square value of the voltage fluctuation, which is obtained experimentally. 


Response of Hot Wire to Three-Dimensional 
Velocity Fluctuations 

Consider the general case when the axis of the hot wire is arbitrarily 
oriented with respect to the mean velocity direction. To anticipate such 
applications in this investigation, it is considered that the mean veloc- 
ity V s is in the 0,z plane, making an angle 0 with the z-direction. 

Let the direction of V s be denoted by subscript s and the direction 
normal to V g by subscript n, the s,n plane being coincident to the 
0,z plane (fig. 37(a)). The stem or holder of the hot wire is always along 
the direction of r. The orientation of the hot wire itself is determined 
by the two angles a and 0 (fig* 57(b))* The angle a is in the 
s,n plane and the angle 0 is between the wire and its stem or the radial 

direction. 
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Equations (Al) or (A2) assume that the steady response of the wire 
varies with the sine of the angle between the wire and the mean-velocity 
V s direction; that is, only the component of V s perpendicular to the 
wire is effective in cooling it. (The validity of this assumption is 
subject to verification by calibration.) In this case, the velocity per- 
pendicular to the wire is Vp = V s sin 7 (see fig. 37(t>)). Expressed in 

terms of the angles a and p, it is 


Vn = Vjl 




p 

sin a 


sin 2 0 


(A4) 


Let the unit vectors along the directions s, n, and r be U 0 . 

— > -» b 

U n , and U r , respectively. Then, the unit vector U qA along the direc- 
tion of the wire is 


— > 
U 


oA 


= sin a sin pu + cos pu + cos a sin pu. 


n 


Hence the unit vector along the direction of V can be expressed by 

Jlt 


Vp _ ? - V^UqA 


V, 


V 


= {1 - sin 2 a Sin 2 p U s - sin a sin P cos P u r . sin a cos a sin 2 p ^ 

]fl - sin^a sin^p {l - sin^a sin^p 

(A5) 

where use has been made of the fact that V s = V s Ug. For smn.1 1 velocity 
fluctuations v s , v n , and v r , only their components along the direction 
of Vp have any influence in changing the magnitude of total Vp and 
hence in cooling the wire. The total fluctuation of V p (due to v s , 
v r , and v n ) as a ratio of itself can be obtained from equations (A4) 
and (A5) as 


dV P _ ^s x vr y vn 
V P " v s " v s " v s 


(A6) 


where 
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sin a sin p cos (3 
1 - sin 2 a sin 2 p 


_ sin a cos a sin^p 
1 - sin^a sin 2 P 


Equation (A6) supplies the value for 


dV 


V T 


of equation (A3) . As a 


check, it is seen that for <x=0 or P =■ 0, the wire is sensitive only 
to v , in agreement with the customary considerations. For p = 90°, 

X = 0 and Y = tan a, and for a = 90°, X = tan p and Y = 0. Both 
facts are well known in hot-wire anemometry. 

Substituting equation (a 6) into equation (A3), squaring and averaging 
over time, gives 


ZlL- + X 2 + Y 2 ZlL. - 2X - 2Y Z& + 2XY ^ = 


V_ 


». 2 


I.1I 2 -I. 2 

Bo ° 


1 + 2| 

(R - ^ 

ii 

l R o J 

•e b 

1 

_ _E_ 

" Eb 



dE 2 

,2 


E 


(AT) 


It is seen that all six components of the turbulent stress tensor 
appear, three components of intensity and three shear correlations. 


Application 

Since the measuring probes used herein are inserted radially through 
the wall of the tunnel, it is quite simple to make use of the previous 
three-dimensional response in the following way: A hot-wire probe is 

made with one support longer than the other j the wire fastened to them 
then assumes just the position shown in figure 37(b). The angle p is 
fixed once the probe is made. Turning the probe stem about its axis has 
the effect of varying a and therefore X and Y. Since there are six 
turbulence quantities to be determined, the wire voltage in six positions 
will be necessary. Two of these were chosen as the positions given by 
a = 0 and a = -90°. The other four are taken as symmetrical about the 
plane formed by the two cross-component fluctuations. Looking radially 
inwards along the probe stem, the scheme is as follows (fig. 37(c)): 
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The two F positions are for a = 0 and should give the same response, 

namely, v s 2 /v s 2 . At the position A the wire is in the r,s plane, with 

its shorter support at the upstream side. This is necessary to avoid 
the wake of this support hitting the wire. The other positions are then 
as shown in figure 37(c). 

Let the right-hand side of equation (A7) he defined as N; that is, 

2 


N = 



1 + 2 5£j dE 2 

1 “ *§ ' e2 


The following system then results: 


V. 




^7 


v 


V c 




Is_ + X 2 ^+Y 2 ^-+2X + 2Y + 2XY ^ = N, 


v e v r 


v^v r 


Vc 


+ x 2 lEI + y 2 — — 2X Iain - 2 Y IsZa + 2XY ^ = N ; 


'D 


+ x 2 + Y 2 IsL _ 2X + 2Y IS*a - 2XY Ilia = N, 


V 0 v v 


V 0 v r 




V c 


Vc 


V 2 
v s 


V, 


E 




+ tan p 


2o v r 


+ 2 tan 0 


v Q v 


s v r 


= n a 


-s 2 


= % 


NACA TN 4135 


44 


Since the coefficients 
except for their sign, this 


in the first four equations are the same 
system is easily solved. It is found that 


v r v n 



+ N D - N B - N e ) 


v v 
s n 




V 


s 


2 


s( NB 

+ N C 

- n d - 

n e) 

5y(" c 

+ n E 

- % - 

n d ) 

Np 




1 

n a - 

% pi 

'n b + n c - n e - n d ^ 

tan (3 

tan 

p x' 

< 4 / 



1 /n b + Nq + Np + N e 

?\ b 



It can be seen from these equations that the presence of turbulent 
shear stresses will be evidenced by an asymmetry of the values of N 
about the F positions. Considering N as a function of a, N would 
be a minimum at the two F positions, if all shear stresses are zero and 
hence Ng = Nq = N d = N e . The increase of N with a is then the same 

amount on both sides of the F positions. If the shear stresses are not 
zero, this is no longer true, and N will not be symmetrical about F. 
The mean voltage or the wire heating current will still be symmetrical 
about F regardless of the turbulent quantities, except for the influ- 
ence of such extraneous factors as possible effects of wire supports, 
uneven plating of the wire, or curvature of the wire. An illustration 
of how heating current I and root -mean-square fluctuating voltage vary 
with angular position is shown in figure 38. These data were taken in 

r 3# £ 

the rotating tunnel, at = 0.35 and — = 17«5* 

b - a M 
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Calibration 

In the usual application of hot-wire anemometry a wire calibration 
is necessary to determine I Q 2 , the extrapolated value of I 2 at zero 

velocity , since this quantity enters in the calculation of N. All other 
quantities that appear in the expression for N are measured directly or 
preset . 

In the present case, it is necessary to know X and Y as well as 
N to calculate the turbulent quantities. The terms X and Y are func- 
tions of a and p, with a determined by the angle through which the 
probe has been turned about its axis and easily measured with a protractor 
arrangement. This leaves 0 and I Q 2 to be found by calibration. This 
is done in the following way. 

The wire is placed in a calibrating airstream, whose velocity can be 
varied, with its axis perpendicular to the velocity (position F). A 
series of points is obtained in the usual way by varying the velocity and 
measuring the wire current. The wire is then turned to position A, which 
implies that the velocity is decreased by cos p if the cosine response 
is valid. A further series of points is then taken. The angle p is 
determined as that angle which gives a straight line for all the points 

on a plot of I 2 against cos p”. An example of this procedure is 

shown in figure 39. It should be noted that the ability to make all such 
points fall on a straight line is an indirect verification of the cosine 
response to angle. With p known, X and Y can then be calculated 
for this particular probe for any value of a. 

It can be seen from figure 38 that a considerable variation of heating 
current is produced by the turning of the wire. Now if only the component 
of velocity perpendicular to the wire axis is effective in cooling, as is 
indicated by figure 39; it is immaterial whether a calibration curve is 
generated by variations in velocity or variations in wire orientation. If 
P is known, the perpendicular velocity component, called V p , may be cal- 
culated for any value of a. Then the variation of I 2 with a can be 
translated to a variation of I 2 with ^Vp, which should be linear. If 

the data of figure 38, for instance, are taken and plotted in this way, a 
straight line does indeed result. It is thus possible, once p is deter- 
mined by a separate calibration as described, to calculate I Q 2 without 
further calibrations from the data obtained during an actual experiment. 

p 

An expression for I Q as a function of heating current for the 

various positions can be derived. The equations for the calibration curve 
are 


46 


NACA TN 4135 


I 


F 


2 



+ I 


2 

o 






2 2 \l/4 

sin a sin $) 





sin^p 


l/4 



The subscripts refer to the various positions of the wire as pre- 
viously described. The slope of the calibration curve is denoted by S. 
From these equations, it is found that 


I 


2 

o 




t 2 

Ib,c,d,e 
sin 2 a sin 2 p) 


174 



Since differences of measured quantities are involved, the accuracy 
in determining I Q 2 from these expressions is not too good unless the 
values of I 2 are very accurate. In practice it is recommended to plot 
a calibration curve from a number of points taken at various values of 
a and to determine I Q 2 from this curve. 

A further check was made for the validity of the cosine law in the 
following way. 

If P = 90°, that is, with a probe whose wire supports are of equal 
length, and if the response to angle a is not cos a but (cos a) 1 !, 
then 

2 / — 2 

i P = S rs + x o 

2 „ r— ri/2 _ 2 

i bi _ S rs cos a i + x o 

2 r )/2 2 

I B2 = S / V s cos a 2 + Z o 
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2 2 

where I-g-^ and. Ig 2 correspond to angles ap and ag- 

From these expressions, 

ip 2 - I B i 2 _ 1 - (cos ai )^ 2 _ c 
I F 2 " I B 2 2 1 - (cos ap ) 11 / 2 

or 


l°g e 

G - 1 

c / cos a l\ n ^ 


\cos 0 * 2 / 

log e (cos a-jJ 


can be derived. 

n "I 

If r] is unity, a plot of ~ against cos a p on log- 

fcOSOn 

G - If i 

| cos (X 2 

arithmic coordinates should yield a straight line of slope 0.5« From a 
number of measurements of 1 ^ as a function of a for a wire with 
p = 90 ° the data shown in figure 40 were plotted. Agreement with the 
line of slope 0.5 seems good. 


Change of Coordinates From r,s,n to r,0,z System 

Once all six components of the Reynolds stress tensor in the 
r,s,n system are obtained, it is a relatively simple matter to find 
these components in the cylindrical r, 0 ,z system on account of the 
relations 

v 0 = v s sin 0 + v n cos 0 
v z = v s cos 0 - v n sin 0 

Hence, 
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V 9 = cos20 + sin 2 0 —— rf + 2 sin 0 cos 0 

V s 2 V V v s 2 



sin^0 -£5 + cos 2 0 — ^ - 2 sin 0 


cos 


v s v n 



v 0 v 


z 



sin 0 



is!) + (=os20 - .ln^^ 
s / s 


V r v 0 


= cos 


0 ^ + 


sin 0 


v r v z 


-sin 0 


V r V n ^ 

■ ■ + cos (ft 


v r v s 



Measurement of Derivatives of Velocity Fluctuations 


For the measurement of the viscous terms in the turbulent energy 
equation it was necessary to measure not only the fluctuating voltage 
signals from the hot wire but also their time derivatives. Using the 
Taylor hypothesis for interchanging time derivatives and space derivatives 
equation (A 7 ) becomes for the wire signal passing through a differentiator 



dv bv bv dv bv dv 

- 2 Y ^ ^ + 2 XY — ^ 
os os OS OS OS OS 



= N' 


where e = dE is the voltage fluctuation 
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In isotropic turbulence the last three terms on the left are zero. 

A careful check based on the symmetry of the signal N' with angle a 
shows that in the flow of this investigation, temns of that type were 
either zero or too small to be measured both with and without rotation. 
Assuming that these terms are zero, the previous procedure is then con- 
siderably simplified, since there are now only three unknowns and there- 
fore measurements at only three positions of the wire are necessary. In 
particular, the above equation becomes 



The transformation from r,s,n coordinates to r,9,z coordinates 
involves derivatives with respect to both s and n but the derivatives 
with respect to n were not measured. To circumvent this difficulty, it 
is assumed here that the isotropic relations 



hold for any one component, though this is not necessarily true for dif- 
ferent components; that is, 



Furthermore, the following isotropic relations are also assumed: 


dv s dv s _ dv n dv n _ ^s ^n _ ^s ^n 
ds dn ds dn ds ds dn dn 


bvs dv n 
Bs ckl 


dy s dvn 

c)n &s 


1 

4 



With these assumptions, it is found that 


(lr) 2 - 
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It is seen that because of the assumption made above the isotropic 
relations still exist for any component in the new coordinates 



These relationships can only be verified by two-wire correlation measure- 
ments, which were not made in the present work. 

The increase of wire signal with increasing deviation of the wire 
from a position perpendicular to the mean velocity offers an interesting 
possibility to check roughly for isotropy in any turbulence field. All 
shear stresses then vanish and the three intensity components are all the 
same . 


Equation (A7) becomes 
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Since 


v s 2 


= Np 


the following relation is obtained: 


1 + X 2 + Y 2 = -JL 

N-p 


Expanding X and Y in terms of a and (3 gives 


1 + 


-.2„ o-i v,2p ^ cog 2 a + cos 2p^ 


N 


(l - sin 2 a sin 2 p) 2 N F 


If 3 = 90°, 


(A8) 


1 + tan^a = X (A9) 

Following a similar procedure for the differentiated signal gives 


1 + 


2 sin 2 a sin 2 p (cos 2 a + cos 2 3) _ IT 
(l - sin2a sin 2 p) 2 


n f * 


(A10) 


if 3 = 90°, 


1 + 2 tan 2 a = (All) 

Np 

N N 1 

In any one of these cases, the ratio — or may be plotted 

Np Np' 

as a function of a, if 3 is known. The actually measured ratios as 
a function of a would indicate isotropy if they follow the relations (A8) 
to (All). As an example, some data taken in the rotating wind tunnel are 
shown in figure 4l, together with equation (A8) calculated for 3 = 62.2° 
associated with that particular wire. The nonisotropy of the flow is 
quite apparent. For comparison, equation (A9) (for 3 = 90°) is also 
shown. It is seen that for 3 < 90° the increase of N with a is less 
than for 3 = 90 °J the reason is physically quite obvious since for small 
values of 3 the wire will remain nearly perpendicular to the mean veloc- 
ity direction for any value of a. 
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Accuracy of. Measurement 

As was previously shown, the determination of the intensities and 
shear stresses involves sums and differences of the measured quantity N 
corresponding to several angular positions of the wire. 

Due to inherent inaccuracies involved in the method of measurement, 
the percentage of accuracy in turbulence level in the mean flow direction 
is believed to be no better than ±5 percent. This uncertainty is reflected 
in the calculations of shear and cross components of intensity; further- 
more, it will depend on a, since errors become critical when the signals 
to be subtracted are very nearly equal. From the response of a wire to a 
turning angle, such as is shown in figure 4l, the following estimate of 
accuracy is made for various values of a for each of the quantities to 
be calculated: 


a, 

deg 

Av r /V 

a, 

deg 

Av n /V 

CL, 

deg 

^sn 

AR rs 

ARnr 

v r /V ' 
percent 

v n /V ' 
percent 

25 

±25 



20 

± 0.08 

± 0.22 


35 

±15 

30 

±22 

30 

±.08 

±.15 

± 0.29 

45 

±10 

40 

±14 

40 

±.08 

±.12 

±.15 

50 

±9 

45 

±11 

45 

±.08 

±.11 

±. l4 

6o 

±7 

60 

±7 

60 

±. o8 

±.09 

±.11 


The table presents the largest possible errors, including inaccuracies 
in N, as well as errors in the determination of the angle. The poor 
accuracy is mainly due to the fact that the change of N with changing 
a is less for p < 90° than for p = 90°. Here it is necessary to use 
P < 90 ° in order to find the shear correlation between the two cross- 
component fluctuations. A penalty is then paid in accuracy of measurement. 
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APPENDIX B 

MOTION OF A PERTURBED ELEMENT IN A ROTATING FLUID 

By Hsuan Yeh 

Consider a main flow field consisting of a solid-body rotation, 
that is, Vg = cur. At a certain instant, say t = 0, imagine a fluid 
element at radius r Q possessing a perturbation velocity in the form 
of an excess tangential velocity v 0 (hence the total velocity of this 
particle is osr 0 + v 0 at t = 0) . The tangential and radial perturbing 
velocities (i.e., velocities relative to the mean rotation) at any sub- 
sequent time will be denoted by v^ and v r , respectively. 

If it is assumed that the pressure fluctuation and viscosity exert 
no influence on the subsequent motion of the fluid element and that the 
angular momentum of the element is thereby conserved, the total tangen- 
tial velocity of this particle at any subsequent radius r will be 
(o>r 0 + v 0 )r 0 /r. The centrifugal acceleration is therefore 

(oor 0 + v 0 ) 2 r 0 2 /r3. However, the pressure gradient of the main flow is 

such that a centripetal force of oAr is imposed at any radius r. 

Hence the net radial acceleration of the fluid particle at any radius r 
is 


dv r _ P r o + v 0 ) 2 r n 2 o 
dt ~ r 3 " “ ' 


It is also possible to derive the same equation based on the fact 
that the total radial acceleration is, in this case, the sum of the 
Coriolis acceleration and the centrifugal acceleration, both due to the 
relative (i.e., perturbed) tangential velocity v 0 . (The centrifugal 
acceleration due to the main flow Vq is balanced by the pressure 
gradient . ) Hence , 


dv 

dt 


r _ 


= 2ldvq + 



But 


= (asp + v Q )r n _ 


oxr 


(B2) 


This again leads to equation (Bl) . 
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Using the familiar substitution 

dv r = d / v r 2 \ 
dt dr\ 2 / 

equation (Bl) can be integrated to become 

2 = .KlLIo)!£o! , a) 2 r 2 + c 

r r2 

The integration constant C can be evaluated by the condition that 
v r = 0 at t = 0 and r = r Q . After a slight rearrangement, the 

result is 



Substituting dr/dt for v r , further integration yields the following 
relation between r and t: 


r 




g(2 + g) 
2 


(i - cos aot) 


(B4) 


y 

where g = — — is a dimensionless ratio representing the strength of 
o>r 0 

initial disturbance. From this the radial as well as the perturbation 
tangential velocities can be obtained as functions of time as follows: 


v r 

(Dr 0 



aor 


o 


g(2 + g) 
2 



g(2 + g) 
2 


g 


g(2 + g) 

2 



g(2 + g) 
2 


sin 2cot 


(l - cos aot) 


(l - cos 2ti>t) 


(l - COS 2CDt) 


(B5) 


(b6) 


These functions are plotted in figure 42 for g = 0.01, 0.05, and- 0.10. 
It is noted that their values do not follow a simple proportion to the 
initial disturbance g. 
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From this simple analysis the following remarks can be made: 

(1) Although the analysis was made for a purely tangential pertur- 
bation at t = 0, the relative velocity subsequently goes through all 

Y 

possible combinations of — . Hence the same analysis covers all pos- 

v e 

sible kinds of initial perturbations, the only difference being a shift 
of the origin t = 0. 

(2) No matter what the initial disturbance is, the subsequent 
(relative) motion of a fluid particle is oscillatory, with a circular 
frequency of 2^>. The maximum radial distance of travel through a com- 
plete oscillation is gr 0 . 


( 5 ) The maximum value of v r occurs at 



which 


point v 0 is zero. Furthermore, ( v r ) ^ = (vq)^ 


max 


( 4 ) The time -averaged value of the product v r v 0 zero * 
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Figure 1.- View of impeller. 
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Figure 2.- Assembly drawing of wind tunnel. 
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Figure 3*- Detail of traversing section. 
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Figure 4.- View of tunnel, looking downstream. 
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L- 57-4344 

Figure 5*- View of tunnel, looking upstream. 
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Figure 6.- Axial velocities without rotation. 
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Figure 7*- Axial velocities with rotation. 


66 


MCA TN 1|-135 



r-a 

b-a 


Figure 8.- Peripheral velocities. 
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Figure 9*- Boundary- layer thicknesses. 
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Figure 11.- Axial velocities with and without rotation. ^ = 17.5* 
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Figure 12.- Peripheral velocity with 4-mesh screens. = 17 *5j 
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Figure 13 •- Turbulence intensities without rotation. — = 17-5. 
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Fig\ire 14.- Shear correlation coefficients without rotation. 
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Figure 15.- Turbulence intensities with rotation. = I7.5. 
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Figure l6.- Shear correlation coefficients with rotation. — = 17 -5* 
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Figure 17 •- Turbulence intensities with rotation. — = 17 -5 • 
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Figure l8.- Shear correlation coefficients with rotation 


NACA TN 4135 


77 




v r v © 

V?" 


1.0 X 10 3 


.8 

.6 

.4 

.2 

0 

-.2 
-.4 
-.6 
-.8 
- 1.0 












n 0 




0 0 


O 

0 

0 



O 

0 

0 




























.2 


.4 


.6 


.8 


1.0 


r- a 
FTET 


Figure 19*- Shear stress distribution with rotation. — = I 7 . 5 . 
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Figure 20.- Effect of rotation on 


NACA TN 4135 


79 


V Z 


.07 


.06 


.05 


.04 


.03 


.02 


.01 












■ 

■ 




o-T 

■ 

□ 

■ 

□ 

■ 

■ 

□ 

□ 

■ 

□ 

■ 

□ n 

□ 


□ 


□ N 

0 ROTA1 
OTATION 

ION 



■ R 





00 .2 .4 .6 

r-a 


.8 1.0 




Figure 21.- Effect of rotation on ^ = 17-5 ■ 
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Figure 23.- Integral scale in mean flow direction 
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Figure 24.- Microscales with rotation. 
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Figure 25 • - Influence of rotation on microscales 
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Figure 27 .- Turbulence decay with rotation. r " - = 0 . 5 . 
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Figure 28.- One -dimensional spectra of turbulence intensities at 

r - a 
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Figure 29.- One -dimensional spectra of v r 2 . 
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Figure 30.- Dimensionless energy balance on v z with rotation. 

= 17*5; PTC denotes the pressure-velocity and triple-velocity 
term; DVW denotes the term combining dissipation and viscous work. 
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Figure 31«- Dimensionless energy balance on Vq^ with rotation. 

= 17 * 5 ; PTC denotes the pressure-velocity and triple-velocity 
M 

term; DVW denotes the term combining dissipation and viscous work. 
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Figure 32.- Dimensionless energy balance on v r 2 with rotation. 

|l = 17«5; PTC denotes the pressure-velocity and triple-velocity 
termj DVW denotes the term combining dissipation and viscous work. 
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Figure 53 •- Dimensionless energy balance on q2, with rotation. 

= 17 * 5 ; PTC denotes the pressure-velocity and triple-velocity 
term; DVW denotes the term combining dissipation and viscous work. 
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Figure 36 .- Turbulent shear stress in a plane channel, from reference 14. 
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(a) Coordinate system. 
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£ direction of mean velocity 
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(b) Generalized position of wire. 
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(c) Hot-wire measuring positions. 

Figure 37«- Sketch illustrating quantities necessary to specify hot- 
wire position and orientation. 
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Figure 38 • - Variation of wire signal with angular position 
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Figure 39*- Typical wire calibration. 
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Figure 40.- Determination of exponent t). G 
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Figure 4l.- Variation of* wire signal with wire turning angle. 
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Figure 1+2.- Variation of tangential and radial perturbing velocities 
with time in solid-body rotation of mean flow. Solid lines indicate 
v r ^/oir 0 j dashed lines indicate Vg^cor 0 . 
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